Next-generation Sequencing Data Research Articles

Interaction of sound-audition traits between eared insects and arthropodophagous bats: using a DNA approach to assess diet

Abstract Arthropod–bat interactions are often considered as a base model for studying factors underlying predator–prey coevolutionary processes. Bats developed ultrasonic echolocation to hunt, and in response some arthropods developed defense mechanisms such as ultrasonic hearing, allowing them to elude bat predators. The present study analyzes the feeding patterns of bats, focusing on sonic-auditory sensory mechanisms in predator–prey interactions. Next-generation DNA sequence data from fecal samples were used to analyze the diet of 17 bat species from Mexico. Arthropod prey taxa were classified according to their auditory traits, and echolocation data were recompiled from literature review. We: (i) classified arthropod families according to their hearing ability; (ii) estimated arthropod taxon richness and proportion in the diet of each bat species; and (iii) used multidimensional scaling, principal component analysis, and regression to analyze prey consumption patterns in relation to their auditory traits and in relation to echolocation characteristics of bats. Finally, we analyzed the relationship between foraging time and auditory characteristics of prey. Families with hearing organs correspond to the orders Lepidoptera and Orthoptera. We registered 20 families of Lepidoptera and 5 of Orthoptera—7 and 3 with hearing organs, respectively. Of these orders, families lacking ears were recorded in the diet of a few bat species. Our results support the allotonic frequency hypothesis predicting a difference in emission frequency intervals between predator and prey. However, we found that the consumption of earless moths is less frequent and is related to diurnal and twilight activity—hence, their consumption is limited to bat species foraging early. Results indicate bats feed on arthropod prey successfully despite the ultrasonic hearing ability of the prey. These results may be due to counteradaptations that allow maintenance of an asymmetric “arms race” between bats and eared insects that favors the predator.

Abstract PO2-08-05: ERBB2 mutated breast cancer detected by next generation sequencing: integrated analysis based on real-world data and The Cancer Genome Atlas

Abstract Background ERBB2 mutation has been recognized as an actionable target for pan-HER tyrosine kinase inhibitors in the era of precision medicine for solid tumors. In our real-world cohort utilizing next-generation sequencing (NGS), we successfully identified cases of ERBB2-mutated breast cancer. Furthermore, we investigated the clinical and pathological characteristics of these cases by integrating our data with The Cancer Genome Atlas (TCGA). Methods Our real-world cohort with NGS included 423 breast cancer patients treated at Gangnam Severance Hospital in Seoul, Republic of Korea, between April 2017 and January 2023. Among these patients, 13 were confirmed to have ERBB2 mutation through NGS analysis. We also incorporated data from TCGA database, which included 818 cases with ERBB2 sequencing. We compared the clinical characteristics (age, sex, stage, and ER, PR, and HER2 status) between the two groups. Furthermore, we examined the frequency of two of the most common mutations, PIK3CA and TP53, in both groups. Additionally, we investigated the distribution of PAM50 subtypes in the TCGA cohort. Results This study analyzed a total of 1,244 cases with NGS data. Among these cases, 31 (2.5%) exhibited ERBB2 mutation, with 13 patients identified in our dataset and 18 in the TCGA data. Among the 31 cases with ERBB2 mutation, 21 (67.7%) had HER2-negative disease. No significant differences were observed in other clinical characteristics based on ERBB2 status. Within the group of 31 ERBB2-mutated cases, TP53 mutations were found in 12 cases (38.7%), while PIK3CA mutations were found in 12 cases (38.7%). In the TCGA dataset, 4 cases (22.2%) with ERBB2 mutation were classified as HER2-E subtype, and the rate of HER2-E subtype was significantly higher in the ERBB2-mutated group compared to the ERBB2-wild group (22.2% vs. 7.6%, p=0.024). Conversely, 14 cases (77.8%) with ERBB2 mutation exhibited a non-HER2 enriched subtype. Conclusions Consistent with previous knowledge, ERBB2 mutations can be detected in HER2-negative breast cancer, despite their rarity. Furthermore, a considerable proportion of ERBB2-mutated breast cancer may depend on various oncogenic signaling pathways other than the HER2 pathway. Further studies are warranted to optimize therapeutic strategies for ERBB2-mutated breast cancer. Citation Format: Junghyun Kim, Yoonwon Kook, Seung Ho Baek, Min Ji Kim, Sohyun Moon, Seung Eun Lee, Jee Hung Kim, Soong June Bae, Joon Jeong, Sung Gwe Ahn. ERBB2 mutated breast cancer detected by next generation sequencing: integrated analysis based on real-world data and The Cancer Genome Atlas [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO2-08-05.

KRAS Allelic Variants in Biliary Tract Cancers.

Biliary tract cancers (BTCs) contain several actionable molecular alterations, including FGFR2, IDH1, ERBB2 (formerly HER2), and KRAS. KRAS allelic variants are found in 20% to 30% of BTCs, and multiple KRAS inhibitors are currently under clinical investigation. To describe the genomic landscape, co-sequence variations, immunophenotype, genomic ancestry, and survival outcomes of KRAS-mutated BTCs and to calculate the median overall survival (mOS) for the most common allelic variants. This retrospective, multicenter, pooled cohort study obtained clinical and next-generation sequencing data from multiple databases between January 1, 2017, and December 31, 2022. These databases included Princess Margaret Cancer Centre, MD Anderson Cancer Center, Foundation Medicine, American Association for Cancer Research Project GENIE, and cBioPortal for Cancer Genomics. The cohort comprised patients with BTCs who underwent genomic testing. The main outcome was mOS, defined as date of diagnosis to date of death, which was measured in months. A total of 7457 patients (n = 3773 males [50.6%]; mean [SD] age, 63 [5] years) with BTCs and genomic testing were included. Of these patients, 5813 had clinical outcome data available, in whom 1000 KRAS-mutated BTCs were identified. KRAS allelic variants were highly prevalent in perihilar cholangiocarcinoma (28.6%) and extrahepatic cholangiocarcinoma (36.1%). Thirty-six KRAS allelic variants were identified, and the prevalence rates in descending order were G12D (41%), G12V (23%), and Q61H (8%). The variant G12D had the highest mOS of 25.1 (95% CI, 22.0-33.0) months compared with 22.8 (95% CI, 19.6-31.4) months for Q61H and 17.8 (95% CI, 16.3-23.1) months for G12V variants. The majority of KRAS-mutated BTCs (98.9%) were not microsatellite instability-high and had low tumor mutational burden (ranging from a median [IQR] of 1.2 (1.2-2.5) to a mean [SD] of 3.3 [1.3]). Immune profiling through RNA sequencing of KRAS and NRAS-mutated samples showed a pattern toward a more immune-inflamed microenvironment with higher M1 macrophage activation (0.16 vs 0.12; P = .047) and interferon-γ expression compared with wild-type tumors. The G12D variant remained the most common KRAS allelic variant in all patient ancestries. Patients with admixed American ancestry had the highest proportion of G12D variant (45.0%). This cohort study found that KRAS allelic variants were relatively common and may be potentially actionable genomic alterations in patients with BTCs, especially perihilar cholangiocarcinoma and extrahepatic cholangiocarcinoma. The findings add to the growing data on genomic and immune landscapes of KRAS allelic variants in BTCs and are potentially of value to the planning of specific therapies for this heterogeneous patient group.

Acute promyelocytic leukemia with torque teno mini virus (TTMV)::RARA fusion: an approach to screening and diagnosis

Acute promyelocytic leukemia (APL) with variant RARA translocation is linked to over 15 partner genes. Recent publications encompassing six cases have expanded the spectrum of RARA partners to torque teno mini virus (TTMV). This entity is likely under-recognized due to lack of clinician and pathologist familiarity, inability to detect the fusion using routine testing modalities, and informatic challenges in its recognition within next-generation sequencing (NGS) data. We describe a clinicopathologic approach and provide necessary tools to screen and diagnose APL with TTMV::RARA using existing clinical DNA or RNA-based NGS assays, which led to identification of four cases, all without other known cytogenetic/molecular drivers. One was identified prospectively and three retrospectively, including two from custom automated screening of multiple data sets (50 257 cases of hematopoietic malignancy, including 4809 acute myeloid leukemia (AML)/myeloid sarcoma/APL cases). Two cases presented as myeloid sarcoma, including one with multiple relapses after AML-type chemotherapy and hematopoietic stem cell transplant (HSCT). Two cases presented as leukemia, had a poor response to induction chemotherapy, but achieved remission upon re-induction (including all-trans retinoic acid (ATRA) in one case) and subsequent HSCT. Neoplastic cells demonstrated features of APL including frequent azurophilic granules and dim/absent CD34 and HLA-DR expression. RARA rearrangement was not detected by karyotype or FISH. Custom analysis of NGS fusion panel data identified TTMV::RARA rearrangements, and in the prospectively identified case, facilitated monitoring in sequential bone marrow samples. APL with TTMV::RARA is a rare leukemia with a high rate of treatment failure in described cases. The diagnosis should be considered in leukemias with features of APL that lack detectable RARA fusions and other drivers, and may be confirmed by appropriate NGS tests with custom informatics. Incorporation of ATRA may have a role in treatment but requires accurate recognition of the fusion for appropriate classification as APL.

Insights into chloroplast genome evolution in Rutaceae through population genomics

Chloroplast genomes, pivotal for understanding plant evolution, remain unexplored in Rutaceae, a family with key perennial crops like citrus. Leveraging next-generation sequencing data from 509 Rutaceae accessions across 15 species, we conducted a de novo assembly of 343 chloroplast genomes, unveiling a chloroplast variation map highlighting the heterogeneous evolution rates across genome regions. Notably, differences in chloroplast genome size primarily originate from large single-copy and small single-copy regions. Structural variants predominantly occurred in the single-copy region, with two insertions located at the single-copy and inverted repeat region boundary. Phylogenetic analysis, principal component analysis, and population genetic statistics confirmed the cohesive clustering of different Citrus species, reflecting evolutionary dynamics in Citrus diversification. Furthermore, a close chloroplast genetic affinity was revealed among Atalantia (previously regarded as primitive citrus), Clausena, and Murraya. Zanthoxylum formed a distinct group with heightened genetic diversity. Through expanding our analysis to include 34 published chloroplast genomes, we explored chloroplast gene selection, revealing divergent evolutionary trends in photosynthetic pathways. While Photosystem I and Photosystem II exhibited robust negative selection, indicating stability, the Nicotinamide adenine dinucleotide (NADH) dehydrogenase pathway demonstrated rapid evolution, which was indicative of environmental adaptation. Finally, we discussed the effects of gene length and GC content on chloroplast gene evolution. In conclusion, our study reveals the genetic characterization of chloroplast genomes during Rutaceae diversification, providing insights into the evolutionary history of this family.

Virus Quasispecies Rarefaction: Subsampling with or without Replacement?

In quasispecies diversity studies, the comparison of two samples of varying sizes is a common necessity. However, the sensitivity of certain diversity indices to sample size variations poses a challenge. To address this issue, rarefaction emerges as a crucial tool, serving to normalize and create fairly comparable samples. This study emphasizes the imperative nature of sample size normalization in quasispecies diversity studies using next-generation sequencing (NGS) data. We present a thorough examination of resampling schemes using various simple hypothetical cases of quasispecies showing different quasispecies structures in the sense of haplotype genomic composition, offering a comprehensive understanding of their implications in general cases. Despite the big numbers implied in this sort of study, often involving coverages exceeding 100,000 reads per sample and amplicon, the rarefaction process for normalization should be performed with repeated resampling without replacement, especially when rare haplotypes constitute a significant fraction of interest. However, it is noteworthy that different diversity indices exhibit distinct sensitivities to sample size. Consequently, some diversity indicators may be compared directly without normalization, or instead may be resampled safely with replacement.

The immune response-related genomic alterations in patients with malignant melanoma.

Immune checkpoint inhibitors (ICIs) significantly improve the survival outcomes of patients with advanced melanoma. However, response varies among from patient to patient and predictive biomarkers are urgently needed. We integrated mutational profiles from next-generation sequencing (NGS) data and clinicopathologic characteristics of melanoma patients to investigate whether tumor genomic profiling contribute to clinical benefit of ICIs treatment. The majority of genes identified with high mutation frequency have all been reported as well-known immunotherapy-related genes. Thirty-five patients (43.2%) had at least 1 BRAF/RAS/NF1 mutation. The other 46 (56.8%) melanomas without BRAF/RAS/NF1 mutation were classified as Triple-WT. We identified mutational signature 6 (known as associated with defective DNA mismatch repair) among cases in this cohort. Compared to patients with PD-L1 expression (TPS < 1%), patients with PD-L1 expression (TPS ≥ 1%) had significantly higher median progression-free survival (mPFS), but no significantly higher durable clinical benefit (DCB) rate. In contrast, FAT1, ATM, BRCA2, LRP1B, and PBRM1 mutations only occurred frequently in patients with DCB, irrespective of PD-L1 expression status. Our study explored molecular signatures of melanoma patients who respond to ICIs treatment and identified a series of mutated genes that might serve as predictive biomarker for ICIs responses in melanoma.

The importance of heteroresistance and efflux pumps in bedaquiline-resistant Mycobacterium tuberculosis isolates from Iran

BackgroundTuberculosis (TB) continues to pose a threat to communities worldwide and remains a significant public health issue in several countries. We assessed the role of heteroresistance and efflux pumps in bedaquiline (BDQ)-resistant Mycobacterium tuberculosis isolates.MethodsNineteen clinical isolates were included in the study, of which fifteen isolates were classified as MDR or XDR, while four isolates were fully susceptible. To evaluate BDQ heteroresistance, the Microplate Alamar Blue Assay (MABA) method was employed. For screening mixed infections, MIRU-VNTR was performed on clinical isolates. Mutations in the atpE and Rv0678 genes were determined based on next-generation sequencing data. Additionally, real-time PCR was applied to assess the expression of efflux pump genes in the absence and presence of verapamil (VP).ResultsAll 15 drug-resistant isolates displayed resistance to BDQ. Among the 19 total isolates, 21.05% (4/19) exhibited a heteroresistance pattern to BDQ. None of the isolates carried a mutation of the atpE and Rv0678 genes associated with BDQ resistance. Regarding the MIRU-VNTR analysis, most isolates (94.73%) showed the Beijing genotype. Fifteen (78.9%) isolates showed a significant reduction in BDQ MIC after VP treatment. The efflux pump genes of Rv0676c, Rv1258c, Rv1410c, Rv1634, Rv1819, Rv2459, Rv2846, and Rv3065 were overexpressed in the presence of BDQ.ConclusionsOur results clearly demonstrated the crucial role of heteroresistance and efflux pumps in BDQ resistance. Additionally, we established a direct link between the Rv0676c gene and BDQ resistance. The inclusion of VP significantly reduced the MIC of BDQ in both drug-susceptible and drug-resistant clinical isolates.

INSaFLU-TELEVIR: an open web-based bioinformatics suite for viral metagenomic detection and routine genomic surveillance

BackgroundImplementation of clinical metagenomics and pathogen genomic surveillance can be particularly challenging due to the lack of bioinformatics tools and/or expertise. In order to face this challenge, we have previously developed INSaFLU, a free web-based bioinformatics platform for virus next-generation sequencing data analysis. Here, we considerably expanded its genomic surveillance component and developed a new module (TELEVIR) for metagenomic virus identification.ResultsThe routine genomic surveillance component was strengthened with new workflows and functionalities, including (i) a reference-based genome assembly pipeline for Oxford Nanopore technologies (ONT) data; (ii) automated SARS-CoV-2 lineage classification; (iii) Nextclade analysis; (iv) Nextstrain phylogeographic and temporal analysis (SARS-CoV-2, human and avian influenza, monkeypox, respiratory syncytial virus (RSV A/B), as well as a “generic” build for other viruses); and (v) algn2pheno for screening mutations of interest. Both INSaFLU pipelines for reference-based consensus generation (Illumina and ONT) were benchmarked against commonly used command line bioinformatics workflows for SARS-CoV-2, and an INSaFLU snakemake version was released. In parallel, a new module (TELEVIR) for virus detection was developed, after extensive benchmarking of state-of-the-art metagenomics software and following up-to-date recommendations and practices in the field. TELEVIR allows running complex workflows, covering several combinations of steps (e.g., with/without viral enrichment or host depletion), classification software (e.g., Kaiju, Kraken2, Centrifuge, FastViromeExplorer), and databases (RefSeq viral genome, Virosaurus, etc.), while culminating in user- and diagnosis-oriented reports. Finally, to potentiate real-time virus detection during ONT runs, we developed findONTime, a tool aimed at reducing costs and the time between sample reception and diagnosis.ConclusionsThe accessibility, versatility, and functionality of INSaFLU-TELEVIR are expected to supply public and animal health laboratories and researchers with a user-oriented and pan-viral bioinformatics framework that promotes a strengthened and timely viral metagenomic detection and routine genomics surveillance. INSaFLU-TELEVIR is compatible with Illumina, Ion Torrent, and ONT data and is freely available at https://insaflu.insa.pt/ (online tool) and https://github.com/INSaFLU (code).

Wide Real-Life Data Support Reduced Sensitivity of Antigen Tests for Omicron SARS-CoV-2 Infections

With the continuous spread of new SARS-CoV-2 variants of concern (VOCs), the monitoring of diagnostic test performances is mandatory. We evaluated the changes in antigen diagnostic tests’ (ADTs) accuracy along the Delta to Omicron VOCs transition, exploring the N protein mutations possibly affecting ADT sensitivity and assessing the best sampling site for the diagnosis of Omicron infections. In total, 5175 subjects were enrolled from 1 October 2021 to 15 July 2022. The inclusion criteria were SARS-CoV-2 ADT combined with a same-day RT-PCR swab test. For the sampling site analysis, 61 patients were prospectively recruited during the Omicron period for nasal and oral swab analyses by RT-PCR. Next-Generation Sequencing data were obtained to evaluate the different sublineages. Using RT-PCR as a reference, 387 subjects resulted in becoming infected and the overall sensitivity of the ADT decreased from 63% in the Delta period to 33% in the Omicron period. This decrease was highly statistically significant (p &lt; 0.001), and no decrease in viral load was detected at the RNA level. The nasal site presented a significantly higher viral load than the oral site during the Omicron wave. The reduced detection rate of Omicron infections by ADT should be considered in the global testing strategy to preserve accurate diagnoses across the changing SARS-CoV-2 variants.

Screening and Diagnosis of Rare Thalassemia Variants.

Rare thalassemia subtypes are often undiagnosed because conventional testing methods can only identify 23 common types of α- and β-thalassemia. To assess a comprehensive approach for the screening and diagnosis of rare thalassemia. The study cohort included 72 individuals with suspected rare thalassemia variants. Screening was conducted by next-generation sequencing (NGS) combined with third-generation sequencing (TGS) and chromosomal microarray analysis (CMA)/copy number variation sequencing. Of the 72 individuals with suspected rare thalassemia, 49 had rare α- or β-gene variants. NGS combined with gap polymerase chain reaction detected a total of 42 cases, resulting in a positive detection rate of 58.3%. Additionally, 4 α-globin genetic deletions were identified by TGS, which increased the variant detection rate by 5.6%. Two samples with a microdeletion of chromosome 16 or 11 were detected by CMA, which increased the detection rate by 2.8%. For one sample, reanalysis of the NGS and TGS data confirmed the presence of the β41-42/βN and βN/βN mosaic. The HBB:c.315 + 2delT mutation was initially reported in Guangdong Province, China. Two HBB gene mutations (HBB:c.315 + 5G>C and HBB:c.295G>A) and 4 rare HBA gene deletions (-11.1, -α27.6, -α2.4, and -α21.9) were initially identified in the Zhonshan region. The hematologic phenotypes of all rare cases in this study were clarified. Rare thalassemia variants are more common than previously thought. Despite advancements in TGS, there is still no foolproof method for detection of all types of thalassemia. Thus, a comprehensive approach is necessary for accurate screening and diagnosis of rare thalassemia variants.

59 Comparison of quality improvement programs relevant in next-generation sequence data analysis

59 Comparison of quality improvement programs relevant in next-generation sequence data analysis

Abstract P19: Comprehensive Benchmark of Somatic Mutation Callers on Cell-Free DNA Next-Generation Sequencing Data in Cancer Liquid Biopsy

Abstract Circulating tumour DNA (ctDNA) has emerged as a rich source of biomarkers in precision oncology. However, ctDNA sequence data is often noisier than that derived from standard tissue biopsies, posing significant challenges for accurate cancer mutation detection. Besides, there is a lack of standardised and unbiased benchmark datasets to evaluate and compare performance of variant callers on ctDNA. To address this issue, we present a novel benchmarking approach and comprehensive benchmark of cancer somatic mutation callers on cell-free DNA (cfDNA) Next Generation Sequencing (NGS) data. To create a large-scale benchmark dataset, we developed a new methodology using real colorectal cancer patients plasma samples exhibiting either ultra-high and an ultra-low levels of ctDNA. This enabled the generation of dilution series of deep Whole Genome Sequencing (WGS, 150x) and the creation of high-quality consensus ground truth calls. We evaluated the performance of 9 SNV callers and 7 short INDEL callers frequently used in academic research, including 2 callers that were specifically designed for cfDNA. We further validated our findings at higher depth of coverage using deep Whole Exome Sequencing (WES, &amp;gt;2,000x) data. We further investigated each caller's features to develop practical guidelines for users depending on their application setting: unguided analysis focusing on overall accuracy vs. guided analysis focusing on sensitivity. Citation Format: Hanaé Carrié, Ngak Leng Sim, Pui Mun Wong, Anna Gan, Yi Ting Lau, Polly Poon, Saranya Thangaraju, Iain Tan, Yoon Sim Yap, Limsoon Wong, Anders Skanderup. Comprehensive Benchmark of Somatic Mutation Callers on Cell-Free DNA Next-Generation Sequencing Data in Cancer Liquid Biopsy [abstract]. In: Proceedings of Frontiers in Cancer Science; 2023 Nov 6-8; Singapore. Philadelphia (PA): AACR; Cancer Res 2024;84(8_Suppl):Abstract nr P19.

Artificial intelligence for Next generation sequencing data analysis

In the rapidly evolving field of genomics, our capacity to decipher genetic data encoded in DNA has been transformed by Next Generation Sequencing (NGS) technologies. These advanced technologies produce an enormous volume of data, posing substantial challenges in extracting meaningful biological insights. Artificial intelligence (AI) algorithms offer distinctive possibilities to unravel the biological information embedded in such extensive and intricate datasets. This review offers a synopsis of AI classifications and algorithms, elucidating how these techniques can be employed on sequencing data. Subsequently, a selection of the most typical, promising, or illustrative applications of AI on NGS data to tackle unresolved technical or biological issues are showcased.

TP53-specific mutations serve as a potential biomarker for homologous recombination deficiency in breast cancer: a clinical next-generation sequencing study

Abstract Background TP53 mutations and homologous recombination deficiency (HRD) occur frequently in breast cancer. However, the characteristics of TP53 pathogenic mutations in breast cancer patients with/without HRD are not clear. Methods Both tumor and paired blood DNA from 119 breast cancer patients were performed clinical next-generation sequencing (NGS) by a 520-gene panel. Mutations, tumor mutation burden (TMB), and genomic HRD scores were assessed from NGS data. Results All TP53 pathogenic mutations in patients were the somatic origin, which was associated with the protein expression of estrogen receptor and progestogen receptor. Compared to patients without TP53 pathologic mutations, patients with TP53 pathologic mutations had higher levels of HRD scores and different genomic alterations. The frequency of TP53 pathologic mutation was higher in the HRD-High group (HRD score≥42) relative to that in the HRD-Low group (HRD score&amp;lt;42). TP53 has different mutational characteristics between the HRD-Low and HRD-High groups. TP53-specific mutation subgroups had diverse genomic features and TMB. Notably, TP53 pathogenic mutations predicted the HRD status of breast cancer patients with an AUC of 0.61. TP53-specific mutations, namely HRD-Low mutation, HRD-High mutation, and HRD common mutation, predicted the HRD status of breast cancer patients with AUC values of 0.32, 0.72, and 0.58, respectively. Interestingly, TP53 HRD-High mutation and HRD common mutation combinations showed the highest AUC values (0.80) in predicting HRD status. Conclusions TP53-specific mutation combinations predict the HRD status of patients, indicating that TP53 pathogenic mutations could serve as a potential biomarker for PARP inhibitors in breast cancer patients.

Towards accurate indel calling for oncopanel sequencing through an international pipeline competition at precisionFDA

Accurately calling indels with next-generation sequencing (NGS) data is critical for clinical application. The precisionFDA team collaborated with the U.S. Food and Drug Administration’s (FDA’s) National Center for Toxicological Research (NCTR) and successfully completed the NCTR Indel Calling from Oncopanel Sequencing Data Challenge, to evaluate the performance of indel calling pipelines. Top performers were selected based on precision, recall, and F1-score. The performance of many other pipelines was close to the top performers, which produced a top cluster of performers. The performance was significantly higher in high confidence regions and coding regions, and significantly lower in low complexity regions. Oncopanel capture and other issues may have occurred that affected the recall rate. Indels with higher variant allele frequency (VAF) may generally be called with higher confidence. Many of the indel calling pipelines had good performance. Some of them performed generally well across all three oncopanels, while others were better for a specific oncopanel. The performance of indel calling can further be improved by restricting the calls within high confidence intervals (HCIs) and coding regions, and by excluding low complexity regions (LCR) regions. Certain VAF cut-offs could be applied according to the applications.

Abstract LB412: BRAFp.G469E hotspot mutation promotes nasopharyngeal carcinoma tumorigenesis with potential druggability

Abstract Nasopharyngeal Carcinoma (NPC) is a distinct type of Epstein-Barr Virus (EBV)-associated head &amp; neck cancer. It is most prevalent in Southeast Asia with a high incidence rate of ~15/100,000. While mutations of the MAPK pathways are found in nearly one-fifth of head &amp; neck squamous cell carcinoma (HNSCC) in Western countries (TCGA-HNSCC Firehose cohort, N=512), and some of which are druggable in HNSCC setting (e.g. MAPK1 and HRAS mutations), we sought to determine the cumulative mutational rate of this potentially druggable pathway in 164 micro-dissected NPC cases using recently published next-generation sequencing data (98 whole-exome and 63 whole-genome-sequenced cases). An overall rate of 9.32% cases of NPC harbored MAPK pathway mutations (16 mutations in 15 tumors). Four genes are recurrently mutated: 4 NRAS (all Q61R/L hotspot mutations), 3 HRAS (G12S, Q61R, S106L), 2 BRAF (R726H, Y647*), and 2 SHC4 (Q519R, A552T) mutations. Among which, both patients bearing BRAF mutations died of NPC, suggestive of the biological importance of BRAF mutations in NPC progression. Prompted by an independent clinical report on an advanced NPC patient with somatic co-mutations of BRAFp.G469E and KRASp.G12D being sensitive to the B-raf/C-raf inhibitor sorafenib, we first examined the biological functions and potential druggability of this particular BRAFp.G469E mutation in NPC, as it represents a hotspot mutation site across 20 cancer types (BRAFp.G469E/A/R/V/* found in 137/70655 pan-cancer samples, cBioPortal). Based on the NPC clinical report, we hypothesized that this BRAFp.G469E hotspot mutation could drive NPC tumorigenesis and might be potentially sensitive to the B-raf inhibitor sorafenib. By ectopic expression of BRAFp.G469E mutation, BRAF-wildtype (wt) and EGFP-control gene into a recently developed Epstein-Barr virus-positive nasopharyngeal carcinoma [EBV(+)NPC] cell line, NPC43, we determined the biological functions of this mutation in NPC. We showed that NPC43-BRAFp.G469E-expressing cells significantly drove a 4-fold increase in 3D tumoroid growth vs. BRAF-wt and EGFP-control cells on day 7 (flat-bottom low attachment plates). Further, wound healing assay demonstrated a significantly higher migration rate of NPC43-BRAFp.G469E cells achieving ~80% wound healing (at 96h) vs. 30% and 10% healing in BRAF-wt and EGFP-control cells. Lastly, our preliminary drug study showed that NPC43-BRAFp.G469E cells displayed higher sensitivity to low dose sorafenib (33.3nM and 100nM) than BRAF-wt and EGFP-control cells, achieving ~20% and 30% growth inhibition at 72h, vs. no growth inhibition in control cells. Thus, BRAFp.G469E hotspot mutation can drive NPC tumorigenesis, and may confer sensitivity to raf inhibitors such as sorafenib in NPC. Potential druggability of other BRAF mutations in NPC should be investigated. Acknowledgements: The study was supported by Research Impact Fund (R4015-19F), Research Grants Council, University Grants Committee, Hong Kong SAR during 6/29/20-9/7/21. Citation Format: Chun Ho Law, Zoey Wing Yee Liu, Franky Leung Chan, Tin Yu Samuel Law, Vivian W.Y. Lui. BRAFp.G469E hotspot mutation promotes nasopharyngeal carcinoma tumorigenesis with potential druggability [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 2 (Late-Breaking, Clinical Trial, and Invited Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(7_Suppl):Abstract nr LB412.

Novel genetic association of the Furin gene polymorphism rs1981458 with COVID-19 severity among Indian populations

SARS CoV-2, the causative agent for the ongoing COVID-19 pandemic, it enters the host cell by activating the ACE2 receptor with the help of two proteasesi.e., Furin and TMPRSS2. Therefore, variations in these genes may account for differential susceptibility and severity between populations. Previous studies have shown that the role of ACE2 and TMPRSS2 gene variants in understanding COVID-19 susceptibility among Indian populations. Nevertheless, a knowledge gap exists concerning the COVID-19 susceptibility of Furin gene variants among diverse South Asian ethnic groups. Investigating the role of Furin gene variants and their global phylogeographic structure is essential to comprehensively understanding COVID-19 susceptibility in these populations. We have used 450 samples from diverse Indian states and performed linear regression to analyse the Furin gene variant's with COVID-19 Case Fatality Rate (CFR) that could be epidemiologically associated with disease severity outcomes. Associated genetic variants were further evaluated for their expression and regulatory potential through various Insilco analyses. Additionally, we examined the Furin gene using next-generation sequencing (NGS) data from 393 diverse global samples, with a particular emphasis on South Asia, to investigate its Phylogeographic structure among diverse world populations. We found a significant positive association for the SNP rs1981458 with COVID-19 CFR (p < 0.05) among diverse Indian populations at different timelines of the first and second waves. Further, QTL and other regulatory analyses showed various significant associations for positive regulatory roles of rs1981458 and Furin gene, mainly in Immune cells and virus infection process, highlighting their role in host immunity and viral assembly and processing. The Furin protein–protein interaction suggested that COVID-19 may contribute to Pulmonary arterial hypertension via a typical inflammation mechanism. The phylogeographic architecture of the Furin gene demonstrated a closer genetic affinity of South Asia with West Eurasian populations. Therefore, it is worth proposing that for the Furin gene, the COVID-19 susceptibility of South Asians will be more similar to the West Eurasian population. Our previous studies on the ACE2 and TMPRSS2 genes showed genetic affinity of South Asian with East Eurasians and West Eurasians, respectively. Therefore, with the collective information from these three important genes (ACE2, TMPRSS2 and Furin) we modelled COVID-19 susceptibilityof South Asia in between these two major ancestries with an inclination towards West Eurasia. In conclusion, this study, for the first time, concluded the role of rs1981458 in COVID-19 severity among the Indian population and outlined its regulatory potential.This study also highlights that the genetic structure for COVID-19 susceptibilityof South Asia is distinct, however, inclined to the West Eurasian population. We believe this insight may be utilised as a genetic biomarker to identify vulnerable populations, which might be directly relevant for developing policies and allocating resources more effectively during an epidemic.

Identification of pine wood nematode (Bursaphelenchus xylophilus) loading response genes in Japanese pine sawyer (Monochamus alternatus) through comparative genomics and transcriptomics.

Pine wood nematode (PWN; Bursaphelenchus xylophilus) is the causative agent of pine wilt disease (PWD), which is considered the most dangerous biohazard to conifer trees globally. The transmission of PWN relies on insect vectors, particularly the Japanese pine sawyer (JPS; Monochamus alternatus). However, the molecular mechanism underlying PWN-JPS assembly remains largely unknown. Here, we found that both geographical and gender could significantly affect the PCA (PWN carrying amount) of JPS; thus, JPS transcriptomes from diverse locations and genders were explored regard to PWN loading. Due to the shortage of genomes, we developed a full-length reference transcriptome for analyzing next-generation sequencing data. A comparative genomic study was performed, and 11 248 potential PWN-carrying associate genes (β) were nominated in JPS by using the reported genomes of PWN and non-PWN carrier insect species. Then, 151 differentially expressed transcripts (DETs), 28 of them overlapped with β, correlated with the PCA of JPS were nominated by RNA-Seq, and found that fatty acid β-oxidation might be the key factor that affected the PCA of JPS. Furthermore, JPS fatty acid β-oxidation rates were experimentally decreased using the inhibitor Etomoxir, leading to an increased PCA of JPS. Meanwhile, silencing MaCPT1 in JPS by RNA interference led to a decreased fatty acid β-oxidation rate and increased PCA of JPS. In conclusion, MaCPT1 was able to decrease the PWN-JPS assembly formation through the fatty acid β-oxidation of JPS. These results provide new insights for exploring the impact of PWN invasion on JPS. © 2024 Society of Chemical Industry.

Clinicopathologic and Molecular Characteristics of HER2 (ERBB2)-Altered Non–Small Cell Lung Cancer: Implications for Precision Medicine

The heterogeneous relationship between protein expression, amplification, and mutations in human epidermal growth factor receptor 2 (HER2) in non–small cell lung cancer (NSCLC) and the optimal methods for detecting these alterations remain unclear. We aimed to elucidate the clinicopathological and molecular characteristics of HER2-altered NSCLC and investigate practical approaches for identifying patients who might benefit from HER2-targeted therapies. Using next-generation sequencing data from 1680 individuals, we searched for patients with HER2-altered NSCLCs, including amplifications and mutations. Clinicopathological data and tissue slides were reviewed. Immunohistochemistry (IHC) and silver in situ hybridization were performed according to the American Society of Clinical Oncology/College of American Pathologists guidelines. Our analysis identified 89 (5.3%) patients with HER2-altered NSCLCs, comprising 30 (1.8%) with amplification and 59 (3.6%) mutations, and they were compared with 165 control patients. Of the 59 HER2-mutated cases, 52 harbored tyrosine kinase domain (TKD) mutations, primarily HER2 exon 20 insertions. HER2 TKD alterations were associated with younger age, female sex, nonsmoking status, adenocarcinoma with a micropapillary pattern, lung-to-lung metastasis, and poor overall survival. The 33 patients with TKD mutations and 3 with non-TKD point mutations showed incomplete or complete membranous HER2 immunoreactivity (1+ and 2+, 61.07%). Six patients exhibiting amplifications had an IHC score of ≤2+ despite their high copy numbers and concomitantly displayed other actionable EGFR, KRAS, SMARCA4, and other HER2 mutations. These HER2-altered NSCLCs with molecular coalterations showed heterogeneous patterns through HER2 IHC and silver in situ hybridization. Therefore, next-generation sequencing should be used to identify HER2 mutations in patients with NSCLC who present with concomitant alterations. In addition, the above clinicopathological characteristics and HER2 IHC results can be valuable determinants for identifying patients with HER2-altered NSCLC. These insights hold promise for the development of more effective diagnostic and therapeutic strategies for this complex subset of NSCLC patients.